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Abstract 

We present a simple quantum-mechanical derivation of correlation function of Langevin sources in the semiclassical 
Boltzmann-Langevin equation. The specific case of electron-phonon scattering is considered. It is shown that the 
assumption of weak scattering leads to the Poisson nature of the scattering fluxes. 
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1. Introduction 

Nonequilibrium electrical noise in mesoscopic systems 
was always the subject of primary interest for Markus 
Biittiker. The famous Landauer-Biittiker formula for the 
shot noise in quantum-coherent conductors became a cor¬ 
nerstone of modern theory of fluctuations [T]. This for¬ 
mula was successfully applied to the calculations of shot 
noise in different mesoscopic systems with noninteract¬ 
ing electrons ranging from double-barrier resonant tunnel 
diodes [5] to quantum-coherent metallic diffusive wires |3] . 
However this method has difficulty in describing interact¬ 
ing electrons or systems with dephasing. To circumvent 
it, one has to introduce dephasing probes laia, i.e. fic¬ 
titious probes with voltages chosen such that they do not 
affect the electrical current but allow a replacement of 
quantum-coherent electrons in the conductor by electrons 
from reservoirs with a random phase. Yet the properties 
of the dephasing probes have to be somehow related to the 
rate of actual microscopic scattering processes. 

An alternative method for calculating the electrical 
noise in conductors in the limit of a large number of quan¬ 
tum channels is the Boltzmann-Langevin approach pro¬ 
posed by Kogan and Shul’man in 1969 In this ap¬ 
proach, the fluctuations of current and any other observ¬ 
able quantities are expressed in terms of the fluctuations of 
semiclassical distribution function, which obey the Boltz¬ 
mann equation with a Langevin source in the right-hand 
side. This method appeared to be very efficient when cal¬ 
culating the hot-electron noise in diffusive metallic wires [7l 
8], frequency-dependent shot noise in metallic structures 
in the presence of external screening gates B EH], and 
even the noise in hybrid superconductor - normal-metal 
systems at voltages much higher than the Thouless en¬ 
ergy m- More recently, it was extended to the case of 
spin-flip scattering in ferromagnetic spin valves |12j and 
applied to Coulomb drag in clean double-layer systems [13] 

The key point in the Boltzmann - Langevin approach 


is the derivation of the correlation function of Langevin 
sources. Kogan and Shul’man derived it assuming that 
the noise arises due to the randomness of electron scatter¬ 
ing by impurities and phonons. It was also assumed that 
all scattering events are independent, hence the scatter¬ 
ing of electrons between a pair of states at a given space 
point presents a Poisson process, whose spectral density is 
proportional to its average rate. 

Surprisingly, there were few attempts to derive the cor¬ 
relation function of Langevin sources directly from quantum- 
mechanical principles. In paper [I4j . this correlation func¬ 
tion was calculated using a sophisticated extension of Keldysh 
diagrammatic technique, which involved time-ordering on 
a four-branch temporal contour. The current paper presents 
a much simpler quantum-mechanical derivation of this quan¬ 
tity, which does not require a diagrammatic technique. 


2. The general expression 


The standard semiclassical distribution function of elec¬ 
trons rip (rt) presents the statistical average of the number 
of electrons in an element of phase space Ap^ x Aa;^ di¬ 
vided by the number of quantum states in this element 
AN = Ap^ Ax^/, which is centred at point (p, r). 
This implies that the statistical averaging is performed on 
top of the coarse-grained averaging [T3]. Once the distri¬ 
bution function is known, one may easily calculate differ¬ 
ent measurable quantities like charge or current density as 
linear functionals of it. It can be shown in many different 
ways USmeidT] that this distribution function obeys the 
well-known Boltzmann equation 




where Icoi is the collision integral that accounts for the 
electron scattering by impurities and phonons or electron- 
electron scattering. This equation is valid provided that 
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the semiclassical approximation holds. First, the charac¬ 
teristic length of spatial variation of np must be larger than 
the microscopic scale responsible for the scattering h/pc, 
where Pc is the characteristic momentum of an electron. 
Second, both the characteristic time of its variation and 
the inverse scattering rate must be larger than the charac¬ 
teristic time of the collision with impurity or phonon 
where Ec is the characteristic energy of an electron. 

Along with the average distribution function, one may 
also be interested in the correlation function of its fluctua¬ 
tions ((5up^ (fiti) 6np^ (r 2 i 2 )), where the fluctuations 6np^ = 
np^ — Up. are only coarse-grained-averaged, and the sta¬ 
tistical averaging is applied to their product [T5]. This 
quantity immediately gives the correlation functions of dif¬ 
ferent observables. To calculate it, let us take a closer 
look at Eq. Q. Apart from the time derivative, the terms 
in the left-hand side describe the deterministic motion of 
electrons in the phase space due to smooth spatial vari¬ 
ations of the distribution function and electrical poten¬ 
tial. In contrast to this, the collision integral describes 
quantum-mechanical transitions of electrons between the 
states with different momentum, which are assumed to be 
local in space and time. These transitions are random and 
should be considered as the source of noise if the semiclas¬ 
sical description is used. 

As the structure of the Boltzmann equation without 
the drift terms resembles the equation of motion of the 
Brownian particle, one may write the corresponding Lange- 
vin equation for the distribution function. To this end, hp 
should be replaced by rip and a random Langevin source 
5J^*{vt) with zero average should be added [H] to the 
right-hand side of Eq. Q . As the duration of an electron 
collision with an impurity or phonon is much smaller than 
characteristic time of variation of Up, this source may be 
assumed to be delta correlated in time. Similarly, it should 
be delta correlated in space because of the local nature 
of the collisions. The momentum-dependent coefficient of 
these delta functions may be calculated as follows [20] . 
Choose an interval of time At much longer than the colli¬ 
sion time fi/Ec but so short that the distribution function 
cannot significantly change during this period. The direct 
integration of the Boltzmann - Langevin equation over 
time gives the increment of Up 
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This is our basic formula for calculating the correlation 
function of Langevin sources. 


3. Equations of motion for the operators 

To carry out the calculations to the end, we have to 
calculate the ratio in Eq. Q using quantum mechanics. 
Consider the particular case of electron-phonon scatter¬ 
ing in an elementary volume of size Ax much smaller than 
the characteristic length at which the average distribu¬ 
tion function or the electrical potential essentially change 
but much larger than h/pc- This allows us to describe 
the scattering with a locally uniform Hamiltonian. The 
Hamiltonian of the system H = Hq -f H is the sum of the 
noniteracting part 

iLo = ^EpO+ap(5) 

P k 

and the part describing the electron-phonon interaction 

= ^^(hk&k + ICk^^k) (®) 

P,k 

where Vk are the matrix elements of electron-phonon in¬ 
teraction, and flp and 6k are the annihilation operators for 
electrons and phonons, respectively. The time-dependent 
occupation-number operator fip = OpOp obeys the Heisen¬ 
berg equation |2T] 


(7) 

As [hp,Ho] = 0, the time dependence of this operator 
is determined only by the weak electron-phonon coupling 
and may be considered as slow. Therefore Eq. 0 may 
be solved by iterations in F. To this end, we perform a 
unitary transformation of all operators 

A{t, t) = i(t) (8) 


l-At 

Arip = np(r,t-|-At)-np(rt) = / drdJ“‘(r,t-l-r). (2) 

Jo 

Hence the correlation function of two such increments is 
given by a double integral 

pAt i-At 

(Anpi(ri) Anp2(r2)) = / dri / dr2 

^0 Jo 

X {Sj;f{r„t + rOSJ;f{r2,t + r2)). (3) 


which brings Eq. Q to the form 

f^ = i|r«.r),6pl. (9) 

If the time interval At is much shorter than the relaxation 
time of the distribution function due to the collisions, the 
increment of fip may be calculated to the second order in 
F, and then the transformation inverse to Eq. ([^ gives 


The delta function of ti — T 2 eliminates one of the inte¬ 
grations, and the other reduces to a multiplication by At. 
Hence it follows from Eq. (|^ that 


Ahp = np{t + At) — fipit) 

l>At 

dr [V{t,T- At), fipit)] 


I 
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X [V{t,T' - At), [V{t,T'' - At), np{t)]], (10) 


where 


p,k 

X e^("p-"p+"'‘)"a++;,k«P (11) 

and all fermionic and bosonic operators are taken at time t. 
The density matrix of the system is assumed to be factor¬ 
ized into the electron and phonon parts, which are diag¬ 
onal in the same representation as Hq. Therefore upon 



averaging Eq. (10), the first summand vanishes and in the 


Figu re 1: (color online). The scattering fluxes that appear in 
Eq. States Pi and p 2 exchange electrons with each other and 

with the rest of states p'. 


second summand, only diagonal terms are left. The av¬ 
erage products of four fermionic operators are decoupled 
into products of pair averages, e.g. “ 

^p-thk.q fiq(l — rip). As a result, all the arguments of ex¬ 
ponents are proportional to the difference r' — t" , and 


^ pAt nlAl 

(Aripi Arip^) =-— J dn J dra 

X {[V{t,Ti - At), npi(t)] [V{t,T2 - At), fip^it)]). (15) 

The averaging is performed according to the same rules as 
in Eq. (10), and the arguments of complex exponents ap- 


rAt 
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dr"Y: 

(^||E_kP cos[(;iw_k + £p+fik - ep)(T' - T")/h] (1 -h A_k) 

-h \Vkf cos[(;iWk + £p - £p+Rk)(r' - t'')/H] Akj 
X Tip (1 

- I |E_kP cos[(Sa;_k + £p+hk - £p){.T' - T")/h] A_k Therefore one obtains 


pear to be proportional to the difference ti — ra. Much like 
when calculating the collision integral, we replace the dou¬ 
ble integration over ti and ra by a sequential integration 
over Ti — Ta and over ri. The former integration transforms 
the exponents in (151 into delta functions of coefficients of 
Ti — Ta, and the latter reduces to multiplication by At. 


+ |VkP cos[(ha;k + - ep^n\c){r' - r")/h] (1 + Ak)| 

^ ^p+/ik (l 

where Ak = (^k^k) are the occupation numbers of phonon 
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+ J(pi -h tik —> pi] — J(pi —^ Pa) — d^(p 2 Pi): (16) 


states. The double integration of Eq. (12) over r' and r' 


can be replaced by a sequential integration over t' — t" 
and t' . If At is much larger than the inverse character¬ 
istic frequencies in the arguments of the cosine functions, 
the integration over t' — t" gives delta functions of its co¬ 
efficients, and the subsequent integration over r' reduces 
to a multiplication by At. Therefore 


(Aup) 


At 

where 
J(p ^ p 

H I 

X (1 -f A_k) + |bkP S{huj]s^ + £p — £p+nk) Afi 


p + tik)], (13) 
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X Tip (1 np_|_;ik) (11) 

is the scattering flux from state p to state p -I- fik. The 


right-hand side of Eq. (13) is exactly the electron-phonon 
collision integral I^oi that appears in Eq. 0- 

Now we calculate the average product of two incre¬ 
ments Ahpj and Ahp^ with the same t and At. To the 
second order in the interaction potential, 


where the scattering fluxes J are given by Eq. ( |14[ ). The 
substitution of this ratio into Eq. (|^ gives precisely the 
same expression for the correlation function of Langevin 
sources that was originally obtained in Ref. [5]. The full 
expression for it is given in Appendix A by Eq. ( |VI| ). The 
physical interpretation of this expression was discussed in 
the above paper and is illustrated in Fig. [T] The arrows 
show the scattering fluxes between different pairs of local 
states, which are independent of each other. If states pi 
and P 2 coincide, all the outgoing and incoming scattering 
fluxes are correlated with themselves and give contribu¬ 
tions proportional to their average rates. If these states 
are different, the contribution is given only by J(pi —> P 2 ) 
and J(p 2 —>■ Pi)- Our calculations show that the inde¬ 
pendence of different scattering fluxes and their Poisson 
statistics are the direct consequences of the assumption of 
weak interaction of electrons with the scatterers. 

The Boltzmann equation Eq. ([^ is nonlinear because of 
the electron-phonon collision integral (13). However in ac¬ 


tual problems, the fluctuations of the distribution function 
are effectively averaged over macroscopic lengths much 
larger than Aa; and therefore are small. Hence the Boltz¬ 
mann - Langevin equation may be obtained from Eq. ([^ 
by linearizing it with respect to the fluctuation drip and 
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adding the Langevin source to the right-hand side 


+ l^-kP (5(Sa;_k -l- £pi — Epj) -^-k 


d5n^ 

~1h 


dSrif 

dr 


f)Sn 

eE^ = 6h-ph + SJ;^K (17) 


Together with the expression for the correlation function 
of Langevin sources given by Eqs. Q and (16), it presents 
the full set of equations for calculating the noise. 


4. Conclusions 

We have derived the correlation function of Langevin 
sources in the Boltzmann - Langevin equation for the fluc¬ 
tuations of the semiclassical distribution function. It was 
postulated that the sources are delta correlated in space 
and time, but the momentum-dependent part was rigor¬ 
ously calculated using quantum-mechanical equation of 
motion for the operators of occupation numbers. If the 
electron-phonon interaction is assumed to be weak, the 
scattering fluxes between different pairs of electron states 
may be considered as independent Poisson processes, and 
the correlation function of Langevin sources is the sum 
of their correlation functions taken with the appropriate 
signs. This coincides with the conditions for the validity of 
the standard collision integral in the Boltzmann equation. 
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Appendix A. 

The full correlation function of Langevin sources for the 
case of electron-phonon scattering is given by the following 
expression: 


{Sj;fir^,h)Sj;f{r2,h))=6ih-t2)Sir,-r2) 

^ 5^1 *^Pi’P2 1^1 + ^pi (1 -^k) 


+ l^-kp ^(£pi+;ik ~ £pi + ^-k) 

^ ^Pi+^k (l ^pi ) 

^Pi,P2 l^^k| ^(^pi ^pi+/ik k) (1 k) 

^(^pi ^pi+k ^^k) -^k 

^ ^Pi "^pi+^iki) 

f^?ik,p2 —Pi *^(^P2 '^Pi ^^k) (1 “1“ -^k) 

+ l^-kP S{hW-\^ + 6p2 — A/"_k 

X (1 TT’pi ) ?T'P2 

*^k.pi—p2 ^(^pi ^P2 ^^k) (1 -^k) 


X np, (1 


-^P2) 


(A.l) 
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